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ABSTRACT 

We present the formulation of a new infinite family of self-consistent stellar models, 
designed to describe axisymmetric flat galaxies. The corresponding density-potential 
pair is obtained as a superposition of members belonging to the generalized Kalnajs 
family, by imposing the condition that the density can be expressed as a regular 
function of the gravitational potential, in order to derive analytically the corresponding 
equilibrium distribution functions (DF). The resulting models are characterized by 
a well-behaved surface density, as in the case of generalized Kalnajs discs. Then, 
we present a study of the kinematical behavior which reveals, in some particular 
cases, a very satisfactory behavior of the rotational curves (without the assumption 
of a dark matter halo). We also analyze the equatorial orbit's stability and Poincare 
surfaces of section are performed for the 3-dimensional orbits. Finall y, we obt a in the 
corr esponding equili brium DFs, using the approaches introduced by IKalnaisl {l976) 
and lDeionghe 



Key words: stellar dynamics - galaxies: kinematics and dynamics. 



1 INTRODUCTION 

The obtention of density-potential pairs (PDP) correspond- 
ing to idealized thin discs is a problem of great astro- 
physical relevance, motivated by the fact that the main 
part of the mass in many galaxies is concentrated in an 
stel lar flat distribution , usua lly assumed as axisymmet- 
ric l|Binnev fc Tremainel ljl987f) h Once the potential-density 
pair (PDP) is formulated as a model for a galaxy, usually the 
next step is to find the corresponding distribution function 
(DF). This is one of the fundamental quantities in galac- 
tic dynamics specifying the distribution of the stars in the 
phase-space of positions and velocities. Although the DF 
can generally not be measured directly, there are some ob- 
servationally accesible quantities that are closed related to 
the DF: the projected density and the light-of-sight veloc- 
ity, provided by photometric and kinematic observations, 
are examples of DF's moments. Thus, the formulation of 
a PDP with its corresponding equilibrium DFs establish a 
self-consistent stellar model that can be corroborated by as- 
tronomical observations. 

Now then, there is a variety potential-densi t y pair s 
for such flat stellar models, e.g. IWvse fc Mavalll l|l942t) : 
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1 20061 ). In particular, IToomrd (|l963lJ 19641 ) formulated a gen- 



eralized fam ily of models w hose first member is the one in- 
troduced bv lKuzminl(|l956l ). This family represents a set of 
discs of infinite extension, derived by solving the Laplace 
equation in cylindrical coordinates subject to appropriated 
boundary conditions on the discs and at infinity. 

Analogously, iGonzalez fc Reinal (|2006l ) obtained a fam- 
ily of finite thin discs (generalized Kalnajs discs) whose first 
member corresponds pr ecisely to the well-known model de- 
rived by IKalnaisl dl972h. S uch family was derived by using 
the Hunter's method ( 196&t ). which is based in the obtention 
of solutions of Laplace equation in terms of oblate spheroidal 
coordinates, by imposing some appropriate conditions on the 
surface density. So, by requiring that the surface density be- 
haves as a monotonously decreasing function of the radius, 
with a maximum at the center of the disk and vanishing at 
the edge, detailed expressions for the gravitational potential 
and the rotational velocity were obtained as series of elemen- 
tary functions. Also, some two-integral distribution func- 
tions for th e first four members of this family were re cently 
obtained bv lPedraza. Ramos-Caro fc Gonzalezl (|2008). Now 
then, as the generalized Kalnajs models correspond to discs 
of finite extension, they can be considered as more realistic 
descriptions of flat galaxies than the Toomre's family. 

In the present paper we formulate a new infinite set of 
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finite thin discs, obtained by superposing the members of 
the generalized Kalnajs family in such a way that the re- 
sulting density surface can be expressed as a well behaved 
function of the gravitational potential. As it was pointed 
out by some authors, this is a fundamental requirement for 
the searching of equilibrium distribution functions ( DF) de- 
scribing such axisymmet r ic syst e ms (see for exampl e, Fricke 
l| 19621 ) . Hunter fc Quianl j 19931 ) . Ijiang fc Ossipkovl fcOO'ft) ). 
Thus, the new family formulated here has the advantage of 
easily providing the corresponding two-integral DFs. 

Furthermore, the models have two additional advan- 
tages. In one hand, the mass surface density is well-behaved, 
as in the case of generalized Kalnajs discs, having a max- 
imum at the center and vanishing at the edge. Moreover, 
the mass distribution of the higher members of the family is 
more concentrated at the center. On the other hand, the ro- 
tation curves are better behaved than in the Kalnajs discs. 
We found that, in some cases, the circular velocity increases 
from a value of zero at the center of the disc, then reaches 
a maximum at some critical radius and, after that, remains 
approximately constant. As it is known, such behavior has 
been observed in many disclike galaxies. 

Now, apart from the circular velocity, there are two 
important quantities concerning to the interior kinemat- 
ics of the models: the epicyclic and vertical frequen- 
cies, which describe the stability against radial and ver- 
tical perturbations of particles in quasi-circular orbits. 
We found that the models formulated here are radi- 
ally stable whereas vertically unstable, which is a char- 
acteristic inherited from the generalized Ka lnajs family 
l|Ramos-Caro. Lopez-Suspez fc Gonzalez! l|2008l )). When we 
deal with three dimensional orbits, there are also a common 
feature between the new models and the Kalnajs family: 
the phase space structure of disc-crossing orbits, that can 
be viewed through the Poincare surfaces of section, is com- 
posed by shape-of-ring KAM curves and prominent chaotic 
zones. However, there are certain situations in which the 
chaoticity disappears and the 3-dimensional motion of test 
particles is completely regular. In such cases, one can sug- 
gest the existence of a third integral of motion, as in the 
case of Stackel and Kuzmin potentials. 

Finally, in order to formulate the new family as a set of 
self-consistent stellar models, we shall deal with the problem 
of obtaining the corresponding equilibrium DFs. By Jeans's 
theorem, they are functions of the isolating integrals of mo- 
tion that are conserved in each orbit. Some authors have 
shown that, if the density can be written as a function of the 
gravitational potential, it is possible to find such kin d of two- 
integr al DFs (see Eddington (1916); Fricke (1952): iKalnaisI 
(1976); Jiang and Ossipkov (2007 )). In this pap er we shall 
adopt the approach introduced bv lKalnaisI (|l976l ). which fits 
quite well to axisymmetric disc-like systems. Then, start- 
ing from the DFs derived from this method, a second kind 
of DFs is obtai ned by using the formulae introduced by 
iDeionghel (1 19861 ). that takes into account the principle of 
maximum entropy. These DFs describe stellar systems with 
a preferred rotational state. 

Accordingly, the paper is organized as follows. First, at 
section [2] we obtain the potential-density pairs for the new 
family of thin disc models. Then, at section we study the 
motion of test particles around these new galactic models 
and, at section [4] we derive the distribution functions asso- 



ciated with the models. Finally, at section [5] we summarize 
our main results. 



2 A NEW FAMILY OF THIN DISC MODELS 

2.1 The generalized Kalnajs discs 

In this subsection we summarize the prin cipal features of the 
genera lized Kalnajs family, introduced bv lGonzalez fc Reinal 
(2006), an infinite family of axially symmetric finite thin 
discs. The mass surface density of each disc (labeled with 
the positive integer n) is given by 



£„(#) = E 



(n) 



i-5 



-1/2 



(1) 



where M is the total mass, a is the disc radius and Ec™' are 
constants given by 



(2n + l)M 
2^2 ' 



(2) 



Such mass distribution generates an axially symmetric grav- 
itational potential, that can be written as 



3=0 



(3) 



Here, P2j(Tj) and Q2j(C) = * 2j+1 Q2j are the usual Legen- 
dre polynomials and the Legendre functions of the second 
kind respectively, — 1 < r\ < 1 and < £ < 00 are spheroidal 
oblate coordinates, related to the usual cylindrical coordi- 
nates (R, z) through the relations 



(4) 



in such a way that the discs are located at £ = and 77 = 
yl — R 2 /a' 2 . Finally, the Cij are constants given by 



C 2 i = 



MG7r 1/2 (4j + l)(2n + l)! 



a2 2 "+i(2j + l)(n - j)T(n+j + §)*y+i(0) ' 



where G is the gravitati onal constant. 

As it was shown bv lGonzalez fc Reinal l|2006l ). in these 
models the surface density is a monotonously decreasing 
function of the radius, with a maximum at the center of 
the disk and vanishing at the edge, being the mass distribu- 
tion of the higher members of the family more concentrated 
at the center. On the other hand, the corresponding rotation 
curves behave as follows: for n — 1, the circular velocity is 
proportional to the radius, whereas for the other members 
of the family it increases from a value of zero at the center of 
the discs, then reaches a maximum at a critical radius and, 
finally, decreases to a finite value at the edge. Besides, the 
critical radius decreases as the value of n increases. 



2.2 Formulation of the New Family 

Now, we will show that it is possible to formulate a new 
family of stellar models by performing a linear combination 
of generalized Kalnajs discs, in such a way that the new 
surface densities can be written as polynomials of the new 
potentials. As it was quoted, this is a basic requirement for 
the derivation of equilibrium DFs through the Kalnajs for- 
malism sketched above. In particular, we are interested on 
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to derive a simple relation between the relative potential on 
the disc and the surface density. 

At first, note that <l?n(0, «), given by ((3j), can be rewrit- 
ten as 



<i>„(0,r?) = -£X> s ,*7 2s ~ 2r , 

3 = r=0 

where A 3r are constants defined as 

(-l) r (4 s -2r)!C 2s g2 3 (0) 



A 3 



2 2s r\(2s - 2r)\(2s - r)\ ' 



(5) 



(6) 



and the relatio n ([SJl was derived by introducing the identity 
(|Arfkenl \200$) ) 



(-l) r (4s - 2r)\ _^ 2s _ 2r 



2 2s r!(2s - 2r)!(2s - r)! 



(7) 



From ([5]) we note that the maximum value of the gravita- 
tional potential on the nth disc is $ n (0,0). Therefore, we 
define the relative potential on such a disc as 



v n (v) = <M0,0)-<M0,n). 



(8) 



Now, suppose that we can choose a linear combination of 
those ^ n leading to a new relative potential ty m of the form 



(9) 



where B n are constants that can be determined from A nr 
and C n (see section l!0)) . 

The new relative potential $ m is generated by a new 
mass distribution described by a surface density E m that is 
also a linear combination of generalized Kalnajs discs E m . 
That is 



(10) 



From this relation and (5), we can note that E m can be 
rewritten as 

( 2 n-l)/( 2 m) 



t m (R) — ^ B n E 



(n) 



(11) 



Thus, we can see that this new family of discs is character- 
ized by the fact that the surface density can be split as a 
combination of powers of the relative potential. This impor- 
tant fact makes viable the further derivation of two integral 
DFs for the whole family (see section [4} , that can be con- 
sidered as a set of self-consistent galactic models. Now, the 
above statements are only true if we can determinate the 
constants B n , introduced in @. In the next subsection, we 
show a procedure that, by using the orthogonality proper- 
ties of P n , leads to a recurrence relation expressing B n in 
terms of A nr and C„. 



where C n are constants defined by 

m 

(0) - 8 n C , 2 l g 2l (0)P 2l (0). (13) 

i=0 

So, by introducing (I12p into ©, ^ m can also be written in 
terms of Legendre polynomials as 

m n m 

* m = ^^ J BnC , l P 2l (r ? ) = ^D„P 2 „(n), (14) 

n=0 i=0 n = 

where D n are constants to be determined. Here, it is impor- 
tant to note that these D n are such that 



Dm — B m Crni (15a) 

D m -i — B m Cm-\ +-B m -iC m _i, (15b) 

Dm-2 = BmCm-2 + Bra-lC m -2 + -B m _ 2 Cm-2, (15c) 



D m -k — C'm-k B m -i. (15d) 

i=0 

The above equations can be written in a compact way as 

m — n 

D n = Cn B m-i, (16) 

i=0 

which summarizes the previous series of recurrence relations. 

Now, let us come to the problem of calculating D n . 
According to ([9]) and ()14|) . we have 



22 D n P2n(ri) = AmO 



n 



(17) 



and by using the orthogonality properties of Legendre poly- 
nomials, we obtain 



D n — AmO 



4n + 1 



J r] 2rn P2n{v)dv, 



that reduces to 

D n = ^4m0 



7r 1/2 (4n + l)r(2m + l) 



2 2m + 1 r(l + m - n)r(m + n + § ) 



(18) 



in such a way that, by using (116[) . we can determine the 
constants B n - In order to do this, note that equations (|15a[) - 
()15d|) give us recurrence relations: from (|15a|l we obtain B m , 
from (|15b|) we obtain B m -i, and so on. In general, we have 

Z) "' 1, (19a) 



B m 
B m ~l 
Bm-2 



Cm 
D m -1 
Cm—1 
Dm-2 
Cm — 2 



Bn 



Bm Bm — 1 , 



(19b) 
(19c) 



2.3 Calculation of B n 

According to definitions (O and ([3]), the relative potential 
associated to the generalized Kalnajs discs can be written 



= ]Ta 1 p 2 „( 7 ?), 



(12) 



B 



m — k 



D„ 



a 



m — k 



i=0 



for k > 1, 



relations that can be summarized as 

m — n — l 



B n = 



Dn 

Cn 



2~2 Bm-i, for n < m — 1. 



i=0 



(19d) 



(20) 
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Figure 1. We plot the mass surface density T,2, £3, £4 and 
S5 (from bottom to top in the left edge), given by (122 b . where 
the parameters B n are determined by J20P . Since the resulting 
densities are negative in some ranges, we must to correct B\. 



m 


Bi 


B 2 


B 3 


B 4 B 5 


2 


-5/8 


1 






3 


-35/192 


-7/12 


1 




4 


-105/1024 


-21/128 


-9/16 


1 


5 


-1155/16384 


-231/2560 


-99/640 


-11/20 1 



Table 1. Constants B n for the first four models: m = 2, 3, 4, 5 



In table Q] we show the values of B n for the first four 
models, i.e. m — 2,3, 4, 5. Although we have solved the prob- 
lem of finding these constants, there is another inconvenient. 
If one introduces such coefficients in (|10p . the correspond- 
ing surface density is negative for certain ranges of R, as it 
is shown in figure [TJ where we plot the resulting E for the 
first four models. However, this problem can be solved by 
correcting B\, i.e. the coefficient corresponding to the dom- 
inant term in (|10|l . In the next subsection, we show that it 
is possible to define minimum values for Bi, corresponding 
to each model, in such a way that all of they be described 
by positive surface densities. 



m 




B 2 


S3 


Bi B 5 


2 





1 






3 


0.101273 


-7/12 


1 




4 


0.128914 


-21/128 


-9/16 


1 


5 


0.143207 


-231/2560 


-99/640 


-11/20 1 



Table 2. Corrected constants B n for the first four models: m 
2,3,4,5 



here formulated, the surface density can be split as 

-,2\ 1/2 



E m (7?,Bi) = BiE 



(i) 



ir 



2 \ n-l/2 



(22) 



where the coefficients B n , for n > 2, are given by (|20[l , A 
simple way to find Bi m in, is by demanding that E m has a 
minimum (equal to 0) at B\ — Bi m i n and R = R m i n . That 
is, we demand that the following two equations holds: 



dE(_R, -Bimin) 



dR 



Bi 



= 0, 



0. 



(23) 



(24) 



The relation (|23p imposes the condition that the surface 
density has a minimum at Bi = -Bimin and R — -Rmin, while 
through the relation (|24|) we demand that its value at such 
critical point is 0. 

We find numerically the lower bounds of -Bi for the 
first four members of the new family and they are showed 
in table [2] The corresponding surface densities, for different 
values of Bi > -Bimin, are plotted in figure [2] We note that, 
in all of these cases, the mass concentration start in for 
R = a and increases towards the disc bulge. Moreover, such 
concentration increases with m. Finally, the gravitational 
potential corresponding to the new family, derived as a linear 
combination of generalized Kalnajs discs, can be cast as 



2.4 Correction of Bi 

From here on, we consider Bi as an arbitrary parameter 
that can be chosen in such a way that E > 0, in the range 
< R < a. For each model, we expect that -Bi has a lower 
limit but does not has an upper limit. The reason is that, 
according to fTJ), we have 



dE m 
hm , „ 

R^o, dR 



-co for m = 1, 
for m > 2. 



(21) 



This means that the behavior of Ei, for R — > a, differs from 
the remaining density surfaces, characterized by a rate of 
change tending asymptotically to at the disc edge. Thus, it 
is evident that one always can find a minimum value -Bimin, 
such that the product BiminEi is larger than any linear com- 
bination of E2, E3, 

In the particular case corresponding to the new models 



(25) 



where B„ are given by (120 ^ , for n > 2 and Bi is an arbitrary 
parameter with lower bound determined by (123|> - (|24 p - 



3 KINEMATICS OF THE NEW FAMILY 

In this section we study the motion of test particles around 
the galactic models formulated above. Since each <J> m is 
static and axially symmetric, the specific energy E and the 
specific axial angular momentum I are conserved along the 
particle motion. This fact restricts such motion to a three 
dimensional subspace of the (R, z, Vr, V z ) phase space. By 
defining an effective potential $J„ as 



2-R 2 ' 



(26) 
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Figure 2. We plot the mass surface density (a) £2, (b) E3, (c) £4, (d) £5, for different values of the parameter B\. In each case we 
start with B\ = B\ m { n (lower curve), given by table[2] The remaining curves, from bottom to top, corresponds to increments of 0.01 in 
Si, starting from B lmin . 



the motion will be d etermined by 
l|Binnev fc Tremaind i|l987l )) 

R = Vr, 

z = V z , 

Vr 



together with 



dR 

dz 



the equations 

(27a) 
(27b) 

(27c) 
(27d) 

(28) 



which gives the total energy of the particle. 

Relations (|27a[l - (|28[) are the basic equations that deter- 
mine the motion of a particle with specific axial angular mo- 
mentum I and energy E, in cylindrical coordinates. At first, 
we restrict our attention on particles belonging the disc, i.e. 
the interior kinematics, in order to describe rotation curves 
(for circular orbits) and the stability of nearly circular orbits 
(by deriving the epicyclic and vertical frequencies). Then, 
we shall focus on three dimensional motion, in particular, 
the case of disc crossing orbits. In such case, we use z — 
surfaces of section, in order to illustrate the regularity or 
chaoticity characterizing those orbits. 

3.1 Interior Kinematics 

We start by pointing out that the system (|27a[) has equilib- 
rium points at Vr = V z = z = 0, R = R c , where R c must 



satisfy the equation 

(Rc,0) 



OR 



R'i 



+ 



dR 



= 0, 



(29) 



(flc.O) 

that is the condition for a circular orbit in the plane z 



0. 

In other words, the equilibrium points of (|27a|l occur when 
the test particle describes equatorial circular orbits of radius 
R c , specific axial angular momentum given by 



4 = ±4 m 



\ 9RJ (R 



(flc.O) 



and specific energy 



E = K l (Rc,0)- 



(30) 



(31) 



The subscript c in i c indicates that we are dealing with 
circular orbits. 

Now, a feature of special interest is the circular veloc- 
ity V c = IcRc (sometimes denoted by V v ), which can be 
directly compared with astronomical observations. Then, it 
is convenient to express V c as a function of the radius. From 
(|30|) . we can write the magnitude of the circular velocity as 

V c 2 = R ( d$ m /dR) . In figure [3] we plot VJR) corre- 
V / z=0 

sponding to the first four models of the new family of discs. 
For Bi = -Bimin the rotation curves has a maximum and 
then decreases to a constant value, in a s imilar fashion that 
the be havior of generalized Kalnajs discs (|Gonzalez fc Reinal 
(2006)). Moreover, the maximum of the rotation curve is at 
an even smaller radius when the value of m increases. Now, 
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Figure 3. We plot the rotation curves for (a) m = 2, (b) m = 3, (c) m = 4, (d) m = 5, for the same values of parameter B\ as in figure 
[2] The lower curves corresponds to B\ m i n and the upper curves corresponds to larger values of B\ . 



in the case in which B\ is very large, the curve approximates 
to a straight line, as a consequence of the dominance of term 
associated with the usual Kalnajs disc. On the other hand, 
for intermediate values of Bi we note an interesting feature: 
when the rotation curve reaches its maximum, it remains 
nearly constant, as it is the case in many real galaxies. 

In order to study the stability of these trajectories un- 
der small radial and vertical (z-direction) perturbations, we 
analyze the nature of quasi-circular orbits. They are charac- 
terized by an epicycle frequency n and a vertical frequency 
v, given by (Binney & Tremaine Il987f ) 

K = U^w ^ = b^W (32) 

This means that by introducing (130[) in the second deriva- 
tives of we obtain k 2 and v 2 as functions of the radius. 
Values of R such that k 2 > and (or) v 2 > corresponds 
to stable circular orbits under small radial and (or) vertical 
perturbations, respectively. In figures [4] and [5] we show the 
behavior of the epicycle and vertical frequency, respectively, 
for the first four models and using the same values of B\ as 
in figure [2] We note that these models are characterized by 
a prominent range of stability under radial perturbations. 
In particular, for B\ > B\ m i n , there will be radially stable 
orbits with radius in the range < R < a. In contrast, there 
are prominent ranges of vertical instability. Such ranges tend 
to decreases when m and B\ increases. Thus, we can say that 
the stability under vertical perturbations, in quasi-circular 
orbits, improves in models with a large m. 



3.2 Exterior Kinematics: Disc-crossing Orbits 

In this section, we study the behavior of 3-dimensional mo- 
tion, i.e. orbits outside the equatorial plane (except when 
they cross the plane z = 0). As it was mentioned above, the 
motion is determined by (|27a[) and can be described in an 
effective phase space with three dimensions (there are two 
integrals of motion: E and I). An adequate tool to investi- 
gate such orbits is the Poincare surface of section, in order 
to find the chaotic and (or) regular regions characterizing 
the structure of the phase space. 

In particular, we present numerical solutions of (|27a[) - 
(I27d|) for the case of bounded disc-crossing orbits. There are 
certain values of E and £ for which they are confined to 
regions that contain the disc and will cros s back and forth 
through it. As it was showed by Hunter (1993), this fact 
usually gives rise to many chaotic orbits due to the disconti- 
nuity in the z-component of the gravitational field, produc- 
ing a fairly abrupt change in their curvatures. There is an 
important exceptional case of this behavior: the Kuzmin's 
disc, characterized by an integrable potential of the form 
$ = -GM{R 2 + (a+ jz|) 2 ]~ 1/2 , with a > 0. However, the so- 
called Kuzmin-like potentials, characterized by &(e) where 
e = [R 2 + (a + |z|) 2 ] 1,/2 , are non-integrable and present the 
behavior mentioned above. The family of models formulated 
here present a very similar structure and we can expect an 
analogous dynamics. Each potential 4 > m(C) 7 ?) can be cast in 
a Kuzmin-like form if we take into account that, according 
to @, £ = (R+ + R-)/2a and -q = (R+ - i?_)/2ia, where 
R+ = [R 2 + (z + ia) 2 ] 1/2 and R- = [R 2 + (z -ia) 2 ] 1/2 . More- 
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Figure 4. Quadratic epicyclic frequency as a function of R for the cases (a) m = 2, (b) m = 3, (c) m = 4, (d) m = 5. We use the same 
values of parameter B\ as in figure [2] The lower curves corresponds to B\ m i n and the upper curves corresponds to larger values of B\ . 
In all of these cases, we note a prominent range of radial stability. Such range increases with B± . 



over, they are characterized by a z-derivative discontinuity 
in the disc, given by 



d<&rr, 

dz 



= + 



d$ m 
dz 



2nG E m (R). (33) 



Despite the above relation makes the KAM theorem inappli- 
cable, we also found a large variety of regular disc-crossing 
orbits. 

In Figure [6] we show the z = 0-surface of section cor- 
responding to some orbits with E = —1.245 and £ — 0.2 
(these same values are used to perform the next three sur- 
faces of section), corresponding to a test particle moving 
around the model m = 2 with B\ — 0.1. As it was ex- 
pecte d, from our experience with the general i zed K alnaj s 
discs (Ra mos-Caro. Lopez-Suspez fc Gonzalez! ((2008)), this 
plot exhibits a variety of regular and chaotic trajectories. 
There is a regular central region conformed by two kinds of 
KAM curves: the central rings made by box orbits; a set of 
resonant islands chain (made by loop orbits) enclosing the 
rings. Moreover, there are two lateral regular zones of loop 
orbits that, as well as the central region, are enclosed by a 
sea of chaos. In figure [7] we show the effect of increasing 
Bi in model m = 2. Then, for the same initial conditions 
with E — —1.245 and £ = 0.2, the surface of section reveals 
an increasing in the chaotic region along with a distortion 
of the KAM curves in regular zones (for example, now the 
central torus are made only by box orbits). 

One could expect similar chaotic surfaces of section for 
models m = 3, . . ., but what really happens is that for cer- 



tain values of Bi all the orbits are regular. This is the case 
illustrated in figure [8] where the Poincare section, corre- 
sponding to a particle moving around the model m = 3, 
reveals completely regular motion. This is a surprising fact, 
since in this kind of potentials chaos is the rule. If one in- 
creases Bi (for example, from 0.2 to 1), as in the case of 
figure [!3 we find again a prominent chaotic region enclosing 
some regular regions of KAM curves. Thus, the free param- 
eter Bi, apart from determining some relevant features in 
equatorial orbits, plays an essential role in 3-dimensional 
motion and we can conjecture that, for certain values, it is 
possible to have a third integral of motion. 



4 TWO-INTEGRAL DISTRIBUTION 
FUNCTIONS 

The derivation of the DFs, associated to the models con- 
structed above, is particularly simple if we work i n an ad- 
equat e rotating frame, as it was pointed out by iKalnaisI 
l| 19761 ). To begin with, we had seen that for a B\ given by 
the recurrence relation (|20[) (from here on we denote it as 
Bi) the relative potential would be 



'.on 



(34) 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

R R 



Figure 5. Quadratic vertical frequency as a function of R for the cases (a) m = 2, (b) m = 3, (c) m = 4, (d) m = 5. We use the same 
values of parameter B\ as in figure [2] The lower curves corresponds to B\ m in and the upper curves corresponds to larger values of B\ . 
The range of vertical stability is small and increases with B\ . 



1.5 




0.2 0.4 0.6 O.B 1 

R 

Figure 6. Surface of section for some orbits with I = 0.2, E = 
— 1.245, around the model m = 2 with B\ = 0.1. 

On the other hand, if we choose a different B\, the relative 
potential becomes 

m 

* m = £1*1 + V B„*„ = A m0 rf m + (Bi - 50*1, 

^2 (35) 

= A m0 ri 2rn + - B^f^a 2 ?? 2 , 

where fi = [37rGM/(4a 3 )] 1/2 . 

Now, let's work in a a rotating frame with angular ve- 
locity Q. In this frame, the effective potential is defined as 



0.2 0.4 0.6 0.8 1 1.2 1.4 

R 

Figure 7. Surface of section for some orbits with £ = 0.2, E = 
— 1.245, around the model m = 2 with B\ = 0.5. 

<|Binnev fc Tremainej i|l987l )) 

$ e = *-^ 2 i? 2 = $ + ^ 2 a 2 (r7 2 -l), (36) 

and by choosing properly the constants, the relative-effective 
potential will be 

= Ano?? 2 " 1 + \{B X - BO^aV - iffaV . (37) 

Note that if one chooses SI as 

n = ±n V B i ~ B i, ( 38 ) 
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ated fi, given by (|38[) , are 



h 



(2) ,.-3/4 



3T,^'A 



(1) 4-1/4 



4JV4 



c V (3) 4-5/6 

6JV6 



4J3/4 



7V( 2 ) 4" 3 / 6 
24J 3 / 6 



+ 



6J 5 / 6 



7^A-^ 45E^A7„ 5/8 



8JV8 



128J 3 / 8 



1024J 5 /8 



+ ■ 



8J 7 / 8 



(42a) 



(42b) 



(42c) 



Figure 8. Surface of section for some orbits with £ = 0.2, E = 
— 1.245, around the model m = 3 with Bi = 0.2. In this case we 
have only regular orbits 



QE (5) A 

JS — 



9/10 



10J 1 / 10 



77E< 4) A 



-7/10 



200J 3 / 10 



,-:!?f.-Vs,j-.. 



4, 




0.2 0.4 0.6 0.8 



Figure 9. Surface of section for some orbits with £ = 0.2, E = 
— 1.245, around the model m = 3 with Bi = 1. Now, a prominent 
chaotic region appears as a consequence of the overlapping of the 
lateral island chains, showed in the previous figure. 



the term with 77 vanishes and equation (|37p reduces to 

* e ,m = A m0 r, 2m . (39) 

Therefore, the relation between the density and the relative- 
effective potential can be written as 



/ f 

— ■ \ A m o 



(2n-l)/(2m) 



(40) 



Finally, by using the Kalnajs method ( Kal naisl (1 19761 ) ). we 
obtain that the DF corresponding to the m-model is given 
by 

m r,(»)/ 9 ,u-(2»-l)/(2™) 

47rm t-^ 



~ (e + ^L, - in"o J ) 1 " (an " 1)/(am) 



(41) 

The explicit DFs for the first 4 models and the associ- 



99Si 3) 4- 5/10 
1280J 5 / 10 



693Sj 2) A 5 - 3/10 
25600J 7 / 10 



1QJ9/10 



(42d) 



where J = e + f2L 2 
45GMtt 



if2 2 a 2 is the Jacobi's integral and 



A 20 = 



128a 



A 30 = 



A 40 



A 6 



175GMrr 
512a : 



11025GMn 
32768a 1 



43659GMtt 
131072a : 



Q 2 3 



3ttGM(5 + 8Bi) 
32a 3 



3ttGM(35 + 192Bi) 
768a 3 



3ttGA/(105 + 1024^! 
4096a 3 



(43a 



(43b) 



(43c) 



37rGM(1155 + 16384Bi) 
65536a 3 



(43d) 

It is easy to see that the DFs obtained by equation (|41[) . 
for the cases m > 3, could be negative in some regions of the 
phase space corresponding to the physical domain. To avoid 
this inconvenient, it is necessary to impose a stronger con- 
dition for the constants B\, in order to obtain well-defined 
distribution functions. To do this, we formulate the following 
equations (in a similar fashion as in subsection 12. 4|l : 



d/ m (J, B' ln 



dR 



0. 



fm (*/min, ^lmin 



: 0. 



(44) 



(45) 



Relation (|44[1 imposes the condition that the DF has a min- 
imum at B\ — B' lmln and R = i? m in, while through the re- 
lation (I45[) . we demand that its value at such critical point 
vanishes. The numeric solution to these equations give us 
the values shown in table[3] for the models with m — 3,4, 5. 
So, taking this values as a lower limit for B\ , the DFs given 
by (|41[) become positive-defined in the physical domain of 
the phase space. The figure [10] shows the graphics of the DFs 
as functions of the Jacobi's integral, with different values of 
B\. In general, we can observe that the probability is max- 
imum for small values of J, and tends to a constant as J 
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m 


Dl 

Imin 


2 





3 


0.182292 


4 


0.287827 


5 


0.381061 



Table 3. Coefficients B' lmin for the first four models: m = 
2,3,4,5 

increases. Moreover, in the cases m > 3 we can see that for 
values of B\ near to -B lmin , the probability has a minimum 
at J ss J m i n , and it is zero for B 1 — B' lmln at J = J m i n , in 
agreement with equations (|44[) and (|45|l . 

On the other hand, it is convenient to derive a new kind 
of DFs, correspond ing to more pro bable rotational states. 
As it was shown bv lDeionghei (|l986l ). it is possible to obtain 
DFs obeying the maximum entropy principle, through the 
equation 

f m (e,L z )= 2 ^ {£ _±\ (46) 
1 + e aJj > 

where f m + is the even part of (|41fl • The figure [11] shows 
the behavior of the DFs given by (146 |l . In (a) and (b) are 
plotted the contours corresponding to the model m = 2, 
with different values of parameter a. As it can be seen, a 
determines a particular rotational state in the stellar system. 
As a increases, the probability to find a star with positive 
L z increases as well. A similar result can be obtained for 



q < 0, when the probability to find a star with negative L z 
decreases as a decreases, and the corresponding plots would 
be analogous to figure 1111 after a reflection about L z — 0. 
In (c) and (d) are plotted the contours of model m = 3 with 
Bi ~ Bimin f° r the same values of a. The behavior of the 
DFs for the remaining cases is pretty similar to the shown 
in these figures: when Bi ^> Bi m i n , the contours are similar 
to (a) and (b), while if Bi ~ -B lm i n , the contours are similar 
to (c) and (d), in agreement with figure [TOl 



5 CONCLUDING REMARKS 

We have obtained a set of models for axisymmetric flat 
galaxies, by superposing members belonging to the gener- 
alized Kalnajs discs family. The mass distribution of each 
model (labeled through the parameter m = 2,3,...), de- 
scribed by (jlOp . is maximum at the center and vanishes at 
the edge, in concordance with a great variety of galaxies. 
Moreover, the mass density can be expressed as a function of 
the gravitational potential (see equation (|22p ). which makes 
possible to derive, analytically, the equilibrium DFs describ- 
ing the statistical features of the models. 

These models have also interesting features concerning 
with the interior kinematical behavior. On one hand, we 
showed that for some values of Bi, the circular velocity has 
a behavior very similar to that seen in many discoidal galax- 
ies. This is a very relevant fact, which suggests that it is 
not always necessary to introduce the hypothesis of dark 
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Figure 11. Contours of fm, given by l|46| l. for m = 2, B\ = 0.1 and (a) a = 1; (b) a = 10. Moreover, it is plotted the case m = 3 with 
B\ =0.2 and (c) a = 1; (d) a = 10. Larger values of the DF corresponds to lighter zones. 



matter halos (or MOND theories) in order to describe ade- 
quately a variety of rotational curves. On the other hand, the 
analysis of epicyclic and vertical frequencies, associated to 
quasi-circular orbits, reveals that the models are stable un- 
der radial perturbations but unstable under vertical distur- 
bances. With regard to the motion of test particles around 
the models formulated here, we found that the behavior of 
disc-crossing orbits is similar to that seen in the generalized 
Kalnajs family. However, for certain values of the parame- 
ter Bi, the Poincare surface of section reveals that one can 
suggest the existence of a (non analytical) third integral of 
motion. 

On the other hand, we find two kinds of equilibrium 



DFs for the models. Such two-integral DFs can be formu- 
lated, at first, as functionals of the Jacobi' s integral, as i t 
was sketched in the formalism developed by Kal naisl l| 19761 ) . 
This class of DFs essentially describes systems which rota- 
tional state, in average, behaves as a rigid b ody. Then, we 
use the procedure introduced by iDeionghd (119860 ■ obtain- 
ing DFs which represents systems with a mean rotational 
state consistent with the maximum entropy principle and, 
therefore, more probable than the first ones. The statements 
exposed above suggest that the family presented here, can 
be considered as a set of realistic models that describes sat- 
isfactorily a great variety of galaxies. 
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